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Su pplementary Figure Legends 

I. Fig. SI: Phenotypic spectrum of bed RNAi cuticles resembles a published allelic series 

A. Unhatched cuticles of embryos laid by MTD-Gal4/UAS-shRNA-bcd females display 
variable numbers of denticle bands. A representative individual from each class is 
shown. The filzkorper, a tail structure, is indicated with arrow heads. Scale bar 200 um. 

B. An additional UAS-shRNA-bcd (GL01320) line gives a similar phenotype. 

C. The maternal driver mat-tub-Gal4 gives a similar phenotype with both UAS-shRNA-bcd 
lines. 

D. To test viability, we arrayed 200 embryos on an agar plate and counted the number 
hatched after 48 hrs at 25C. For bed RNAi embryos, 0/200 and 0/200 embryos hatched; 
for embryos laid by MTD- Gal4/UAS-shRNA-GFP females 222/239 and 178/200 of 
embryos hatched (TRiP Toolbox stock 182) (Neumuller et al., 2012). 

II. Fig. S2: Mother age contributes strongly to the strength and variability of the bed RNAi 
phenotype while temperature, UAS-shRNA-bcd line, maternal driver, zygotic UAS-shRNA- 
bcd construct copy number, and paternal genotype do not. 

A. The severity of bed phenotype increases as the flies ages. To quantify the effect of 
maternal age, we allowed flies to eclose for 48 hrs and collected cuticles from one cage of 
MTD-Gal4/UAS-shRNA-bcd females every few days and counted denticle bands. Old 
mothers are >10 days old and used to build the gene expression atlas. Counts: Day 2 n = 
137, Day 6 n = 173, Day 8 n = 224, Day 11 n = 86, Day 14 n = 127, Day 16 n = 58. 

B. qPCR indicates >90% knockdown in embryos laid by older mothers. Embryos were 
collected for 2 hrs. As a reference for WT bed levels we averaged deltaCT counts from 
Days 3, 5, 8 and 12 of MTD-Gal4/UAS-shRNA-GFP embryos. We then subtracted this 
average number from the bed RNAi deltaCT counts for Days 3, 5, 8, 10, and 14, 
converted the deltadeltaCT to percent of WT levels. 

C. Paternal genotype does not meaningfully influence the bed RNAi phenotype. MTD-Gal4l 
UAS-shRNA-bcd virgin females were crossed to males homozygous for UAS-shRNA-bcd 
or with males homozygous for a enhancer lacZ reporter (WT). Progeny from the first 
cross will have 1-2 copies of the UAS-shRNA-bcd construct (blue), while progeny from 
the second cross will have 0-1 copies of the UAS-shRNA-bcd construct (red). We could 
not detect a difference between the number of ventral denticle bands visible in each 
populations (p values from Kolmogorov-Smirnov test). Note that the effect of mother 
age is much greater than paternal genotype (compare left and right for Day 1 and Day 6 
samples respectively). Although it is possible that shRNAs against bed are zygotically 
expressed, they do not meaningfully contribute to the phenotype, consistent with the 
purely maternal effect of bed. This cross also shows there is no detectable paternal effect. 
The absence of zygotic or paternal effects allowed us to integrate the reporter data into 
the gene expression atlas. 

D. Two UAS-shRNA-bcd lines yield similar phenotypes. We tested two shRNA lines 
(GL00407 and GL01320) at 25°C and 29°C, taking samples on Day 4 and Day 10. For 
these crosses we used the mat-tub-Gal4 maternal driver; this driver tends to give a more 
consistent phenotype than MTD-Gal4. ANOVA of temperature, U AS -line, and day reveal 
that each has a significant but small effect, less than one segment in all cases. For the 
mat-tub-Gal4 driver, the distribution of phenotypes laid by young mothers approaches the 
steady state distribution seen for old mothers with MTD-Gal4 more quickly. This result 



is consistent with more uniform phenotypes for the mat-tub-Gal4 driver with shRNAs 
against other genes (Staller et al., 2013). The distributions of embryos laid by old 
mothers of both genotypes are comparable. For future work depleting other maternal 
effect genes, we plan to use the mat-tub-Gal4 driver. 

III. Fig. S3: bed RNAi embryos have more cells and altered cell density patterns. 

A. Average cell density maps of WT and bed RNAi embryos. While the physical shape of 
the embryos remains asymmetric, the posterior density pattern is duplicated in the 
anterior of bed RNAi embryos, like some mRNA patterns. Embryos from stage 
5:51-100% (time points 5 and 6 in the gene expression atlases) are shown. 

B. Histogram of cell counts in bed RNAi embryos and WT (transgenic) embryos. 

IV. Fig. S4: The boundaries of the eve stripes move over stage 5 in bed RNAi embryos. 

A. Boundary positions calculated using inflection points of WT and bed RNAi embryos. 
Stripe 4 and 5 appear in a handful of embryos is cohort 4, but not frequently enough to 
reliably quantify boundary position. In each plot, anterior is left, dorsal is top. 

B. The widths of the eve stripes contract in bed RNAi embryos. At T = 5 eve stripes 4-7 are 
approximately 1 .7, .6, 1 .4, and 1 .3 cell widths wider in bed RNAi than WT. At T = 6 eve 
stripes 4-7 are approximately 1 .3, .6, 1 .5, and .3 cell widths wider in bed RNAi than WT. 
Data calculated from the side of the embryo. 

V. Fig. S5: The coefficient of variation of most the gap and eve stripe mRNA pattern widths 
are similar between WT (blue) and bed RNAi (red). The exceptions are Kr and the anterior 
eve stripe, which are more variable in bed RNAi embryos. The ventral region of the ectopic 
anterior hb pattern (DV strips 7-11) is very faint in bed RNAi embryos, and our analysis 
script struggles to reliably find a boundary, so this analysis likely overstates the variability 
in this region. Pattern widths calculated with the inflection point. The ectopic anterior hb 
pattern in bed RNAi is compared to the hb posterior pattern in WT. 

VI. Fig. S6: The locations of each transcription factor combination in each genotype. The 
"sum" of each class is shown at the bottom. 

VII. Fig. S7: Changing the ON/OFF threshold does not meaningfully change the conclusions of 
the combination analysis. 

A. For a range of thresholds and most time points, all combinations present in bed RNAi 
also present in WT. 

B . When a combination is unique to bed RNAi, this is generally because it was not detected 
in WT for that threshold and time point. At high thresholds, overlap between adjacent 
patterns {Kr and kni or kni and gt) are not detected in WT. At T=4, the adjacent hb, hkb 
and til patterns do not overlap enough to be detected in WT at most thresholds, but this 
combination is found at adjacent time points, so it is not a true new combination. T=6 is 
explained in the main text. 

C. Line traces of Kr, hkb, and til at T=6 in WT and bed RNAi. Anterior-posterior position is 
on the x-axis and expression level is on the y-axis. The high levels of background in the 
WT hkb pattern may confound the combination analysis. 



VIII. Fig. S8: Thersholds, stain hapten, and mother age do meaningfully influence the fraction of 
cells that express both eve and ftz. 

A. The hapten (DIG or DNP) of the eve mRNA probe does not bias this analysis./te is 
always in the opposite channel. 

B. A schematic shows the effects of varying the ON/OFF threshold for either eve or ftz. To 
find the threshold we create a histogram of the expression level of each gene separately 
and identify the peak of the OFF cell population (mode). For our normal threshold, we 
add one standard deviation. We varied the threshold to be 0.6, 0.8, 1 and 1 .2 standard 
deviations. 

C. C-F) Regardless of the ON/OFF thershold used, the fraction of cells expressing both eve 
and ftz (double ON cells) decreases over time. The mixed mom eve DNP embryos 
(magenta) were collected from cages that had a wide range of mother ages, before we 
started collecting from aged cages. We included these data because they have a very 
high quality stain. These data indicate that mother ages does not have a meaningful 
effect on this analysis. In Fig. 6 the red and magenta data are combined. In all plots, 
error bars are the standard error of the mean. 

IX. Fig. S9: In bed RNAi the eve2 enhancer reporter has limited expression of stripe 7 in the 
posterior and some duplicated stripe 7 in the anterior. 

A. A maximum intensity projection of the eve2 lacZ stain. 

B . DN A (green) , ftz (blue) , and lacZ (red) . 



Table S 1 : The numbers of individual embryos for each gene at each time point included in the 
bed RNAi gene expression atlas 

Table S2: The standard deviations of each gene in the bed RNAi gene expression atlas. Data 
from the WT atlas is included for comparison (Fowlkes et al., 2008). 

Table S3: ON/OFF thresholds used in the combination analysis in Fig. 5, S7. 

Table S4: Sequences of the enhancer lacZ reporter constructs included in the bed RNAi atlas. 
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Identies and Counts of Gap Gene Combinations Unique to bed RNAi 
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Table S1: Embryos per gene bed RNAi atlas 20140106 



Genes 


1 


2 


3 


4 


5 


6 


cad 


5 


15 


24 


14 


3 


6 


D 


9 


19 


26 


13 


5 


13 


eve 


5 


10 


24 


19 


28 


13 


fkh 


5 


14 


16 


9 


4 


7 


ftz 


38 


67 


81 


24 


14 


25 


gt 


25 


30 


31 


13 


8 


9 


hb 


16 
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13 


17 


10 
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hkb 


16 


11 


20 


14 
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14 
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Kr 
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16 


10 
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10 


hb posterior enhancer 
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20 
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3 


gt posterior enhancer 
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18 
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eve3+7 enhancer 
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10 


13 
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10 


eve4+6 enhancer 
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27 
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16 
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eve5 enhancer 
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eve4+6mini enhancer 
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12 
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3 


run 


8 


20 


22 


16 


9 


5 


tn 


15 


20 


20 


11 


6 


6 


hbP 


7 


10 


14 


7 


6 


9 










Average 


11.4 


18.1 


23.8 


12.8 


7.6 


8.8 


Sum 


217 


343 


453 


243 


144 


167 



Table S2: Atlas Standard Deviations 
20140228 



Gene Name 


WT standard bed RNAi standard 
deviation deviation (DNP) 


V/CIU 


0.165 


0.0664 


p\/p 


0.129 


0.1206 


fkh 


0.068 


0.0268 


ft7 

I LZ. 


0.131 


0.1689 


nt 


0.108 0.0717 


hh 

1 IU 


0.134 


0.0585 


hkh 

1 1 r\kJ 


0.106 


0.0365 


kni 


0.099 


0.0583 


Kr 


0.066 


0.0813 


til 


0.063 


0.0365 


D 




0.1009 


run 


0.1509 


hbProtein 

hb posterior 

enhancer 

gt posterior 

enhancer 

eve3+7 

enhancer 

eve4+6 

enhancer 

eve5 

enhancer 




0.1069 




0.0590 




0.0771 




0.0741 




0.0873 




0.1247 


eve4+6mini 
enhancer 




0.0742 



Table S3A: WT thesholds 





T = 1 


T = 2 


T 1 

T = 3 
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T =5 


T C 
1 =D 


fit 


u.oo 
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U.OH- 
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n 91 


Kr 
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n 9^ 
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u.o^ 
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U.OH- 
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til 


0.29 


0.24 


0.19 


0.18 


0.17 


0.16 


hkb 


0.24 


0.30 


0.29 


0.27 


0.23 


0.21 


hb mRNA 


0.41 


0.28 


0.23 


0.23 


0.28 


0.25 


hb Protein 


0.21 


0.20 


0.22 


0.22 


0.22 


0.16 


Table S3B: bed RNAi thesholds 




T = 1 


T = 2 


T = 3 


T = 4 


T =5 


T=6 


gt 


0.27 


0.30 


0.31 


0.26 


0.22 


0.17 


Kr 


0.24 


0.31 


0.31 


0.32 


0.26 


0.24 


kni 


0.20 


0.23 


0.27 


0.29 


0.25 


0.20 


til 


0.26 


0.28 


0.24 


0.19 


0.15 


0.10 


hkb 


0.15 


0.21 


0.24 


0.26 


0.19 


0.18 


hb mRNA 


0.21 
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0.25 


0.22 


0.23 


hb Protein 


0.08 
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0.14 


0.24 


0.27 


0.35 



Construct 



Table S4: Enhancer list 

Enhancer sequence 



Reference Primer 



evel 



eve46 



eve5 



eve3+7 



eve2 



eve46mini 



hb 

posterior 



gt 

posterior 



eve late 
seven 



AGGCCTAATCACTTCCCTGAAATGCATAATTGTGCCGCGGCTTTTGATACGCTCCTGGCGGAGAGGGAGATGAGGAAAGGATGCACGGGAACCGC 

AGCCAAGTGGCAGTCGAGATTGGCAAATCCGCCAGCGGACAATGCCCAGAGAATGGGCAACAAGTAGCGGCGAATTAGCAATCCTATCATGCTTTT 

ATGGCCGGCCAACTCTTGCCCGCGCATCTCAGTTCATCCGAAGCGGGACCAGGTCCAGGTTCAAGTCGAGGTCCAGTACCCCTGCTATCCCGTCA 

ACCCCTTTAGGGCGATAATCCTTCTAAATGTTTGCATTAATTTCGAGGCGTGGACGGATTAGGGCGTGCTGGCTGGGCGGAACCCGCAGCAGAAAC 

CGCCGAGGACACTGCACCGACTGACCTGCAGCCTACAGATCTCTGATCTTCGATCTCTAATCCTTTCGCATTTGCAACTGACTTCTGCACTGGGTCC 

GCCCCTAATCCTTCCGCCGAGAAGGCGGCAGAGTCGCGAGGTACTGGCCCGGGGTAATGGGATTATCTGCGATTACCCCAGATGATCCGCAGAAA 

GTCAATCTGGTTCAGGGGCTAATTGTCAGCGAAGTCAACTAAATCCAATCCTTTCGCGCCCCCTTCTGTTTATTTGTTTGTTTTCGTTTGTTTTGAGA 

ATTTCTGGCAATTAAGTTGCCCGTTTTGATGCGCGGGGGCGGGTGCATCAAATCCTTTCGGCATACCTGTCCTGCACAAATGCTGAATTCCGCATC 

CCATGGATACCCAGATATTCAGATATCCCAAGGC 



Fujioka et 
al., 1999 



AGGATCCCTGGGCTCTGGGCTCTGGACTATCCGCCGACCCTCCATATCCATGATTTACAATTCTCGTTTTTTTCGCGTTATTTTTTTAGGGGCTTTAA 

TGACCGTCGTAAAGCCGCAGGAGGACCAGGACCAGGACTCTGCTCACATTTCGCGCACTGATTCTAAAAAATGAAATCATTTTTTCTTGAATTTCAC 

GGCGCGCCTCGAGCAGGACTCTTTGTTCTCGGCCAGGCAATTGTCCTTTTTTGCGCTCAGCTCTCAGTTTTTTCGTCCAGCGGGCATTACCTACAC 

GGCGTTTTATGGCGGAGATGATATTCGCCTGGGATCGGTTCCGTTTTTTAGGCCATAAAAATTAGGCGGCATAAAAAAACTGCATTGGAATTCTAGT 

TCTAGTTTCAAGTTTTTAGGTTTCCAGGTTTCTGCCAGCCCGCCTAGATTCGCATTTCGCGGAATTCGGAAGCGGAACAGAATGCCAGAATGGTCAG 

AATCCTGGCTGACCTTGCCTTTTGGCCAGGGGCCGTAAAAAAATTGACTCGCTGCGGTGCGCGGAATATTTTTTAAATCTGACTTTCCAACAATCTC 

TGATCTGGGTTCGAATCGTAAAAAAAAAGCAGAACAAAAAGCGGGCATTTTCGTCGGCAAATGATCTGTTAATGGGCCGGGCTAAAAAACTAAGTCA 

CAAAGTCACAAGGTTGTCCGGTAAATTGACCCGGTTAAGAATGTCTGTCTGTACCGAGAAGGATGCAGGACATTCAGCACTTCAAAGCTCCCACCG 

CTCGAAGGATTCCCCCGAAGATTCAC 



Fujioka et al. 
1999 



sequence 

AGGCCTAATCA 
CTTCCCTG 



GCCTTGGGATAT 
CTGAATATCTGG 



ATATCCCAAGGCCGCAAAGTCAACAAGTCGGCAGCAAATTTCCCTTTGTCCGGCGATGTGTTTTTTTTTTAGCCATAACTCGCTGCATTGTTTGGGC 

CAAGTTTTTCTTCTGCCAAATTGCGGAGATGATGCGGGGATTATGCGCTGATTGCGTGCAATTATGGACATCCTGCGAGGCCCCGAGGAACTTCCT 

GCTAAATCCTTTCATCCGCCTACAGAACCCCTTTGTGTCCCGTTCGCCGGGAGTCCTTGACGGGTCCTTCGACTATTCGCTTACAGCAGCTTGCGT 

AAAATTTCATAACCCTACGAGCGGCTCTTCCGCGGAATCCCTGGCATTATCCTTTTTACCTCTTGCCAATCCGTTGGCTAAAAAACGGCTTCGACTTC 

CGCGTAACTGCTGGACAACAAAGACAAAAAACGGCGAAAGGACGGCGATTTCCAGGTAGCATTGCGAATTCCGTCAAACTAAAGGACCGGTTATAT 

AACGGGTTTATATGGCCAGAATCTCTGCATCTCCACGACCGCCAGAAGCTGCGTAAAACTGCAGGCTCTGTTTTGATTTCTGCAACTTCAGTTAATT 

GCCCGGGATGGCCAGCAATTGCCGGCAATTATAAAACAGCGCAGATGTGACTCAGCTTCCATATCTAACTCTATATCTCATGCCGAAAATCTAGGGT 

GGGGAGCGGAGGGGCGGGGTGCGTGGGTGACTTGCCTGCCAGGGAAAGGGGGCGGGGGTTCAGCGGGTGATAAATGTGCGTGATTTGGAATGA 

ATG CG C ATCG ATTAAAACCG C AGG G C AATC AATTT 



GGATCCTCGAAATCGAGAGCGACCTCGCTGCATTAGAAAACTAGATCAGTTTTTTGTTTTGGCCGACCGATTTTTGTGCCCGGTGCTCTCTTTACGG 

TTTATGGCCGCGTTCCCATTTCCCAGCTTCTTTGTTCCGGGCTCAGAAATCTGTATGGAATTATGGTATATGCAGATTTTTATGGGTCCCGGCGATC 

CGGTTCGCGGAACGGGAGTGTCCTGCCGCGAGAGGTCCTCGCCGGCGATCCTTGTCGCCCGTATTAGGAAAGTAGATCACGTTTTTTGTTCCCAT 

TGTGCGCTTTTTTCGCTGCGCTAGTTTTTTTCCCCGAACCCAGCGAACTGCTCTAATTTTTTAATTCTTCACGGCTTTTCATTGGGCTCCTGGAAAAA 

CGCGGACAAGGTTATAACGCTCTACTTACCTGCAATTGTGGCCATAACTCGCACTGCTCTCGTTTTTAAGATCCGTTTGTTTGTGTTTGTTTGTCCGC 

GATGGCATTCACGTTTTTACGAGCTC 

GGTTACCCGGTACTGCATAACAATGGAACCCGAACCGTAACTGGGACAGATCGAAAAGCTGGCCTGGTTTCTCGCTGTGTGTGCCGTGTTAATCCG 

TTTGCCATCAGCGAGATTATTAGTCAATTGCAGTTGCAGCGTTTCGCTTTCGTCCTCGTTTCACTTTCGAGTTAGACTTTATTGCAGCATCTTGAACA 

ATCGTCGCAGTTTGGTAACACGCTGTGCCATACTTTCATTTAGACGGAATCGAGGGACCCTGGACTATAATCGCACAACGAGACCGGGTTGCGAAG 

TCAGGGCATTCCGCCGATCTAGCCATCGCCATCTTCTGCGGGCGTTTGTTTGTTTGTTTGCTGGGATTAGCCAAGGGCTTGACTTGGAATCCAATC 

CCGATCCCTAGCCCGATCCCAATCCCAATCCCAATCCCTTGTCCTTTTCATTAGAAAGTCATAAAAACACATAATAATGATGTCGAAGGGATTAGGG 

G 



TCGAGCAGGACTCTTTGTTCTCGGCCAGGCAATTGTCCTTTTTTGCGCTCAGCTCTCAGTTTTTTCGTCCAGCGGGCATTACCTACACGGCGTTTTA 
TGGCGGAGATGATATTCGCCTGGGATCGGTTCCGTTTTTTAGGCCATAAAAATTAGGCGGCATAAAAAAACTGCATTGGAATTCTAGTTCTAGTTTCA 
AGTTTTTAGGTTTCCAGGTTTCTGCCAGCCCGCCTAGATTCGCATTTCGCGGAATTCGGAAGCGGAACAGAATGCCAGAATGGTCAGAATCCTGGC 
TGACCTTGCCTTTTGGCCAGGGGCCGTAAAAAAATTGACTCGCTGCGGTGCGCGGAATATTTTTTAAATCTGACTTTCCAACAATCTCTGATCTGGG 
TT 

TAGCACGAAAAACCGAAGGATTAAAAAAGGAAACTAGAGCAGAGGTC^ 

CGCCGTAAAATTGGCTATGCGGCCAAACAATAGTGCGAAGGACGACGGCAGGACGCGCAGGACAATCGTCTGGTGGATTTCCAGTCGACACGCCACG 
AGATTTTATGAAGGCAACTCGCTTTGCATGTTATTCCATAGATTTCGCTTCGGTCCCGGTTTGTTTTGGTCAGGTAAGA 

GCTGGAAAATCGCGAGAAACTTCGAAAGACACACAAAGATACAATATCTATGAGTCTAATGGTCATTAGAGCGGTGCGCTCTACATACAATTGTACCAGC 
CGTCTTGTTTGAAGCCTAAAAAACGTCGCAAAAAACACACTTCCGCGTAAGACATCCCATTTCTGTGGTCCGATCGTAAAATATTTAG I I I I I I ATGACCA 
ACGGTGCGGGCAGGTAGCTGGCTGCCG I I I I I I GTGCGCGACCTCAACCCTTTCACCCATTAAGAAAAAATCGCATCCTGTGAGTGTCCTTGCCCGTTT 
CCCTCGAAACGGCCCACAATTTGTGTGCTTTGCGTTTTCTCCTCTC I I I I I GTTTCCACCTAATGTCGGCGTCATTGTCTTCTTTATGACGCCTCGGTTGT 
TTCI I I I I I ATGGTGTCCTTTGTCCTTTGAGCCTCGTTGCACGGCCAAATCCCTACTTCCTCAACCCTTTGGCGGACGAGAAAGTTGCTAGGAGGAGAAC 
GGGTTAAGCGAAAACTCCATTGCACTTTTTACAAGCCGCGATCTTCTTGGAATTAGTTTTGGTCATTAGGCGAAAGGGT^ 
TGGGTTTACTGAGTGAATTCAATGGGCTAAGGCGAGTAAAGGGTTATACTGTTTTTACAT^ 
CCAGCACTTTTAAAAGCCATATATAACTTTAT^ 

GTCTAAAATTTC I I I I I I I I AAATAA I I I I I AAA ATA I I I I I I I GATAGCATACGAAGTATTTAAAAATGTGAACAGATTAAACACATTAAATTTATAAAAGTAAA 
TACAACAGATTTAGCATAGAAATAAAAATCATTTTAATGTTCCGTC 

TTGGCCGAACATTTGGTGCGATTACATTCGTAATTCGCTGGAAATTAAGGCCACTAAGTCGCCAGCGAAATGAATTCGGACATTGGGCATTG 

TAAAAAGGACTCTAGAGCCCCGACCATTGCAATGGTCCATTGTTGAGCGTCCGAAAGATCTGAAAACCAAACCCAAACCAAATCCCGAGCTTAGGCAAT 

CGGCATTGGGAAATAAGCGCCAAATATTCTACCCCCCACTCCAAAAACGAGCATT 

AATTT ATTACCAGAAACTTACCATCACTTCGAGATGAAAGTGCGGAGGAGAGCGCTGATCTTTGCTTATCCCAAACTAGGGGTTTTTTAGGGGTAAA 

AAGCGGGAATACAGATTCTCAGAAACAAAATGGCAGTCCTAATATGTGATAAGTTGCTTTCAACTTACGTTTACGACAGGGGCCATTCGAAGCTAAG 

GATTCCAATATGCGACTGTTAACCCATTAAGACAAAGGGCCGCGAAAGGAGTTAACGGCAATTCCCCGAAACCAAGCACGAATCCAAATCCAAACC 

ACTCCTCACCCTTTTTTGGACGGGTGGGACGGGTCATATAAGGCAAAACCCCTTTCAGTCAATTAGTCGCTGATTTTCTTTGTCACCTAGCAGCGGA 

CCAATATAAAAATCGCAGCCACAATGGTCGGAGGAGAGAACCCTTTTTTTTTTAAGGACCGCCGGTGTCCGAAATATCAGTTTATGGCTCCTTAAAA 

AACTGCAGCGGTTTTAGGGCCCGCGGACTCGGAATGAGGCCTTTTCGCACGAAGCGTCCATTATGCAATAAAACTTTGGATGTTTTTTAGGCAACA 

GCATAACAGCCTAATACAGCGCATACCGTGGGTGGGCAACCCGTTCGGCCATCAGGTAGATGATTCCGTAAAAAACGCGCCCAAGGTCACAACTC 

AAAGGATTGCACAAATATTGCATCTTTGGTTCGGAGCTCATTATGGCGAAGGAACAGCAGTCCGAGTATATGAATATATGGTTTTATGGCCTGCTGT 

TTTTTACGTTTTTTTTTCGCTGAATCTGGTTTTTACACCGAACTTGGCAACTCTTTGGCTTCGACTTCGGCCGCGACTTTTTATGGCATCTGCTCGGG 

ATCGATTGGGCTGGGCTATACGTATGTCCAGGATCTGGAAGGGCTCGGGTTCGGGCCCGGTAACCGCAGGTAAGGGATCGACTCAATGACGGCG 

ACGTGACGGTCCGAGTTTTAGTTTAATCCACCATTTTTACGGTTGGATTAGCGCACGGATTAGCGGATAAGTTCGCGGTTTTTTACTGATTACCATCG 

ATCGATCCCTCTTTTTATTTGCAATTCACCCAATCCCCCTGAATGGGCAGTTGTAAA 
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1 Supplemental Note 1: Patterning by an exponentially dis- 
tributed morphogen 



The following model is reproduced from Eldar and Barkai 2003 and chapter 8 in Alon 2007. 
It assumes that an exponentially distributed morphogen is the only regulator setting the 
position of a boundary. In the D .melanogaster embryo, maternally deposited bed mRNA 
is translated to form an exponentially distributed protein gradient (Gregor et al. 2007). In 
WT, the anterior boundary of Kr sits at 44% of egg length from the anterior and shifts to 
27% of egg length from the anterior in bed RNAi.The anterior boundary of kni shifts from 
59% egg length in WT to 46% egg length in bed RNAi embryos. There is evidence that Bed 
protein activates Kr and kni, and may set the anterior positions of these boundaries. 




100 



Percent Egg Length 

Figure 1: An exponentially distributed morphogen patterns a boundary. The WT gradient 
is in blue and the perturbed (RNAi) gradient is in red. In this figure, f3\ — 1, /?2 = 0.5, 
and a — 0.05. 

We begin with an exponentially distributed morphogen bcdwT(X). 

bcd WT (X) = f3 ie ~ aX (1) 
bcd RNAi (X) = f3 2 e~ aX (2) 

In this example, the boundary of interest is set at the threshold 7. The RNAi perturbation 
only changes the input level of the morphogen from f3± to /?2- The length scale of the 
gradient, a is unchanged. It remains unclear if the bed gradient forms through passive 
diffusion or through a combination of active transport and diffusion, but for this analysis, 
all that matters is that the shape of the gradient be exponential and that the perturbation 
does not change the length scale. This model further assumes that the gradient will be 



1 



exponential regardless of the total input level, an assumption supported by experimental 
data (Liu et al. 2013). 



X\ is the location of the pattern boundary in WT and X2 is the location of the pattern 
boundary in the perturbed embryo. 

7 = fae-** 1 (3) 

7 = /3 2 e- QX2 (4) 

The shift in boundary position is: 

AX = X\ — X 2 (5) 
To solve for AX we set equations (3) and (4) equal to each other. 

7 = f3 ie ~ aXl = f3 2 e~ aXl (6) 

log(f3 ie - aX ) = log(f3 2 e- aX ) (7) 

logtfi) - aX 1 = log((3 2 ) - aX 2 (8) 

log(h) - log{fo) = aX 1 - aX 2 = a(X 1 - X 2 ) (9) 

MA)-MM =Xl _ X2 (10) 



a 



AX = -log(^) (11) 

Importantly, the shift in the boundary, AX is independent of the initial position, Xi, 
and depends only on the relative change in input level of the morphogen, It follows 
that all boundaries set by this morphogen will shift by the same amount, AX. If this 
simple patterning mechanism was responsible for setting the boundaries of Kr and /cm, the 
boundaries of these genes would both shift by the same amount. If Bed were setting both 
the anterior and posterior boundary of a pattern, that pattern would shift, but not expand. 
In reality, the anterior Kr boundary shifts by 17% egg length and the kni boundary shifts 
by 13% egg length. The fact that each boundary shifts by a different amount implies other 
regulators help set these boundaries (see Discussion). 
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